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Q ' Abstract 

O' 

£NJ ■ Exact characteristic trajectories are specified for the time-propagating Wigner phase-space distribu- 

tion function. They are especially simple — indeed, classical — for the quantized simple harmonic oscilla- 
tor, which serves as the underpinning of the field theoretic Wigner functional formulation introduced. 
Scalar field theory is thus reformulated in terms of distributions in field phase space. Applications to 
. duality transformations in field theory are discussed. 



oo : 

An autonomous formulation of Quantum Mechanics, different from conventional Hilbert space or path 
integral quantization, is based on Wigner's phase-space distribution function (WF), which is a special 
representation of the density matrix In this formulation, known as deformation quantization Q, 
phase-space c- number functions are multiplied through the crucial non-commutative ^-product Q. The 
empowering principle underlying this quantization is its operational isomorphism M to the conventional 
J> , Heisenberg operator algebra of quantum mechanics. 

Here, we employ the ^-unitary evolution operator, a "^-exponential" , to specify the time propagation of 
Wigner phase space distribution functions. The answer is known to be remarkably simple for the harmonic 
oscillator WF, and consists of classical rotation in phase-space for the full quantum system. It thus serves 
as the underpinning of the generalization to field theory we consider, in which the dynamics is specified 
through the evolution of c-number distributions on field phase space. 



Wigner functions are defined by 

f(x,p) = ^jdyp(x- \y) e'^x + |y) . (1) 

Even though they amount to spatial auto-correlation functions of Schrodinger wavefunctions ip, they can 
be determined without reference to such wavefunctions, in a logically autonomous structure. For instance, 
when the wavefunction is an energy (E) eigenfunction, the corresponding WF is time-independent and 
satisfies the two-sided energy ★-genvalue equations || , 

H*f = f*H = Ef, (2) 

where * is the essentially unique associative deformation of ordinary products on phase-space, 

* = e \{dxd p —8 p dx) ^ (3) 
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as denned by Groenewold [[|, and developed in [Q]. In practice, it may be evaluated through translations 

of function arguments, f(x,p) k g(x,p) = f(x + ^ dp, p — y dx) g(x,p), to produce pseudodifferential 
equations. 

These WFs are real. They are bounded by the Schwarz inequality H to — 2/h < f < 2/h. They 
can go negative, and, indeed, they do for all but Gaussian configurations, so they are not probability 
distributions Q. However, upon integration over x or p, they yield marginal probability densities in p and 
x-space, respectively. They can also be shown to be orthonormal 0, |5[. Unlike in Hilbert space quantum 
mechanics, naive superposition of solutions of the above does not hold, because of Baker's || fundamental 
nonlinear projection condition f k f = f /h. 

Time-dependence for WFs was succintly specified by Moyal through the commutator bracket H bearing 
his name, 

d 

ih—f(x,p;t) = H*f(x,p;t) - f(x,p;t)*H. (4) 

This turns out to be the essentially unique associative generalization of the Poisson Bracket |?J , to which 
it reduces as Ti — > 0, yielding Liouville's theorem of classical mechanics, d% f + {/, H} = 0. 



For the evolution of the fundamental phase-space variables x and p, time evolution is simply the 
convective part of Moyal's equation, so the apparent sign is reversed, while the Moyal Bracket actually 
reduces to the Poisson Bracket. That is, the ^-dependence drops out, and these variables, in fact, evolve 
simply by the classical Hamilton's equations of motion, x = d p H, p = —d x H. 



What is the time-evolution of a WF like? This is the first question we address. Relying on the 
isomorphism to operator algebras of indicated, one may solve for the time-trajectories of the WF, which 
turn out to be notably simple. By virtue of the ^-unitary evolution operator, a "^-exponential" 0, 

U*(x, p; t) = el tH/n = 1 + (it/h)H(x, p) + W£t H *H + * H -k H + (5) 

the time-evolved Wigner function is obtainable formally in terms of the Wigner function at t = through 
associative combinatoric operations completely analogous to the conventional formulation of quantum 
mechanics of operators in Hilbert space. Specifically, 

f{x,p;t) = U+ l {x,p;t)-kf(x,p;0)-kU+(x,p;t). (6) 

As mentioned, the dynamical variables evolve classically, 

dx x -k H — H -k x 
dt ih 

and 



d P H , (7) 



<g = p*H--H*p = 
dt in 

Consequently, by associativity, the quantum evolution, 

x (t) = U+tx + U- 1 , (9) 

p{t) = U i ,*p*U~ 1 , (10) 
turns out to flow along classical trajectories. 



What can one say about this formal time-evolution expression? Any WF in phase space, upon Fourier 
transformation resolves to 

f(x,p) = f dadb f(a, b) e iax e ibp . (11) 
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However, note that exponentials of individual functions of x and p are also ★-exponentials of the same 
functions, or ★-versions of these functions, since the ★-product trivializes in the absence of a conjugate 
variable, so that 

e iax e ibp = e mx e ibp. ( 12 ) 



Moreover, this is proportional to a ★-product, since 

ibp ia(x+ih c) p/2) ibp iax ibp —ihab/2 



(13) 



Consequently, any Wigner function can be written as 

f(x,p) = Jdadb f(a,b) e lhah ' 2 e^*e^. (14) 
It follows then, that, by insertion of U± ★ U+ 1 pairs at every ★-multiplication, in general, 

f(x,p;t) = f dadb f(a,b) e ihab/2 e ^ _1 *^ * jbu^+p+u* 

dadb f(a, b) e ihab ' 2 eT^ * . (15) 



Unfortunately, in general, the above steps cannot be simply reversed to yield an integrand of the form 
f(a, b) e *»(-t) e «!P(-t). But, in some fortuitous circumstances, they can, and in this case the evolution of 
the Wigner function reduces to merely backward evolution of its arguments x,p along classical trajectories, 
while its functional form itself remains unchanged: 

f(x,p;t) = f(x(-t),p(-t);0). (16) 



To illustrate this, consider the simple linear harmonic oscillator (taking m = 1, u = 1), 

p 2 + x 2 x — ip x + ip h 
= 2 ^2 * V2 + 2 ' ( ' 

It is easily seen that 

ihx = x ★ H — H ★ x = ihp , ihp = p* H — H -kp = —ihx , (18) 
and thus the canonical variables indeed evolve classically: 

X = x(t) = Ui, ★ x ★ U^ 1 = x cos t + psint, P = p(t) = ★ £>★ U^ 1 = pcost — xsini. (19) 
This checks against the ★-exponential for the SHO, Q, ^tH/h _ ^ 1 exp(2itan(t/2)H/h). 

Now, recall the degenerate case of the Baker-Campbell-Hausdorff combinatoric identity for any two 
operators with constant commutator with respect to any associative multiplication, thus for any phase- 
space functions £ and n under ★-multiplication. If 

£*r7-r7*£ = c, (20) 

then, 

eUel = e^fl 2 . (21) 
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Application of this identity as well as ( |i"3| ) and ( |12| ) yields directly 

iax(-t) ibp(-t) ifoab/2 i{p, cos t+bsin t)x+i(b cos t— a sin t)p 

_ i(a cos t+fe sin t)a: i(6 cos t-a sin t)p ih(a cos t+o sin t) (6 cos t-a sin t) /2 
i(a cos t+6 sin t)x cos t— a sin t)p 

= e* e* 

g*(i cos t+fe sin t)x gi(6cos 4— a sin t)p (22) 

Consequently, 

t) = J dadb /(a, 6) e iax ^ e ibp ^ , (23) 
and hence the reverse convective flow (|l6|) obtains. 



The result for the SHO is the preservation of the functional form of the Wigner distribution function 
along classical phase-space trajectories: 

f(x,p; t) = f(x cos t — psini, pcos t + x sint; 0). (24) 

What this means is that any Wigner distribution rotates uniformly on the phase plane around the origin, 
essentially classically, even though it provides a complete quantum mechanical description. Note how, in 
general, this result is deprived of import, or, at the very least, simplicity, upon integration in x (or p) 
to yield probability densities: the rotation induces shape variations of the oscillating probability density 
profile. Only if, as is the case for coherent states ||lQ|| , a Wigner function configuration has an additional 
axial x — p symmetry around its own center, will it possess an invariant profile upon this rotation, and 
hence a shape-invariant oscillating probability density. 

The result (Mj, of course, is not new. It was clearly recognized by Wigner It follows directly from 
(|) for (0) that 

(d t +pd x -xd p )f{x,p;t) = . (25) 

The characteristics of this partial differential equation correspond to the above uniform rotation in phase 
space, so it is easily seen to be solved by (p4|) . The result was given explicitly in Q and also ||, following 
different derivations. B Lesche fl^| , has also reached this result in an elegant fifth derivation, by noting 
that for quadratic Hamiltonians such as this one, the linear rotation of the dynamical variables (|H]) leaves 
the symplectic quadratic form invariant, and thus the ^-product invariant. That is, the gradients in the 
^-product may also be taken to be with respect to the time-evolved canonical variables (JToj) , X and P; 
hence, after inserting U+-kU+ 1 in the *- functional form of /, the ^-products may be taken to be with respect 
to X and P, and the functional form of / is preserved, (|i~6|). This only holds for quadratic Hamiltonians. 



Dirac's interaction representation may then be based on this property, for a general Hamiltonian con- 
sisting of a basic SHO part, Hq = (p 2 + x 2 )/2, supplemented by an interaction part, 



H = H + Hi. (26) 



Now, upon defining 



w = icf + e?" 10 '*, (27) 

it follows that Moyal's evolution equation reads, 

d 

i/i— w(x,p; t) = 7ii * w(x,p; t) — w(x,p; t) * TLj, (28) 
where Hi = ^Ho/h ^ ^ ^ e ^ itH /h^ Expressing Hj as a ★-function leads to a simplification. 
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In terms of the convective variables (|19|), X, P, Tii(x,p) = Hj(X,P), and w(x,p;t) = f(X,P;t), while 
•k may refer to these convective variables as well. Finally, then, 

ih^ f(X, P; t) = H T (X, P) * f(X, P; t) - f(X, P; t) * H T (X, P). (29) 

In the uniformly rotating frame of the convective variables, the WF time-evolves according to the interaction 
Hamiltonian — while, for vanishing interaction Hamiltonian, f(X,P;t) is constant and yields (24). Below, 



in generalizing to field theory, this provides the basis of the interaction picture of perturbation theory, 
where the basis canonical fields evolve classically as above []. 

To produce Wigner functionals in scalar field theory, one may start from the standard, noncovariant, 
formulation of field theory in Hilbert space, in terms of Schrodinger wave-functionals. 

For a free field Hamiltonian, the energy eigen-functionals are Gaussian in form. For instance, without 
loss of generality, in two dimensions (x is a spatial coordinate, and t = in all fields), the ground state 
functional is 

= exp (--L J dx 4> (x) ^m 2 - V 2 (a:)) . (30) 

Boundary conditions are assumed such that the i/m 2 — V 2 kernel in the exponent is naively self-adjoint. 
"Integrating by parts" one of the y/m 2 — V 2 kernels, functional derivation 5(j> (x) /5<f> (z) = 5 {z — x) then 
leads to 



, 5 2 



8<f) (w) 54> (z) 



= ^m 2 - V 2 (j) (w)^j (^Jm 2 - V 2 z 0(z)^J - ti^jm 2 -V 2 z 5(w - z) (32) 



Note that the divergent zero-point energy density, 



E = ^limJm 2 -V 2 5(w - z) , (33) 
may be handled rigorously using ^-function regularization. 

Leaving this zero-point energy present, leads to the standard energy eigenvalue equation, again through 
integration by parts, 

l -Jdz ^-h 2 j^ + <p(z)(m 2 -Vi)cj ) ( Z ^m=E m- (34) 



A natural adaptation to the corresponding Wigner functional is the following. For a functional measure 
[drj/2-K] = Yl x dr] (x) /2tt, one obtains 



W[<f),ir] 



dr] 



2vr 



h 

2 V 



-i J dx r)(x)ir(x) 



h ' 

+ 2^ 



(35) 



where n (x) is to be understood as the local field variable canonically conjugate to 4>(x). However, in 
this expression, both <p and ir are classical variables, not operator-valued fields, in full analogy to the 
phase-space quantum mechanics already discussed. 



1 Ref discusses a field theoretic interaction representation in phase space, which does not appear coincident with the 
present one. 
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For the Gaussian ground state wavefunctional, this evaluates to 
~dr)~ 



2vr 



exp 



— J dx (d (x) - -r) (x) \ \j m 2 - V 2 X (d> (x) - ^77 ( 



(36) 



1 



n 



x e-^/^^M^expj-^-y dx U(x) + ^r](x)\ y>m 2 -V 2 U( x ) + '^i 1 {x) 



n 



<\\-p j — — / dxcf) (x) \j m 2 — V 2 (x) 



dr) 



2vr 



-i J dx r)(x)ir(x) 



exp \ — — j dxrj (x) \J m 2 — V 2 7] (x) ] ) . 



So 



W[<f>,ir] = M exp (~ J dx M (x) ^m 2 - V| <j> (x)j + ^tt (a?) (j/m 2 - 7r(a?)^J, (37) 

where is a normalization factor. It is the expected collection of harmonic oscillators. 



This Wigner functional is, of course ||, an energy *-genfunctional, also checked directly. For 

H [(j>, tt] = ~ J dx (vr (x) 2 + (z) (m 2 - V 2 ) <p (x)) , 
and the inevitable generalization 

ih f . I 8 5 5 5 



x) 5ir(x) 5tt(x) 5(p(x) 



(38) 



(39) 



it follows that 
H *W = 



dx 

T 

dx 



( m2 " v ') ^TxT + * (x) ( m2 " v ') * (x) " I 



2^5tt (x) 



4 (57T (x) 



4 (50 (x) 



T,^[0,7r] 



£ W [0, vr] . 



(40) 



This is indeed the ground-state Wigner energy-*-genfunctional. The *-genvalue is again the zero-point 
energy, which could have been removed by point-splitting the energy density, as indicated earlier. There 
does not seem to be a simple point-splitting procedure that regularizes the star product as defined above 
and also preserves associativity. 

As in the case of the SHO discussed above, free-field time-evolution for Wigner functionals is also 
effected by Dirac delta functionals whose support lies on the classical field time evolution equations. Fields 
evolve according to the equations, 



— ifidt4> = H -k 4> — <ft* H, —ifidtit = H -kir — it * H. 

For Ho, these equations are the classical evolution equations for free fields, 

dt4> (x, t) = tt (x, t) , dtTT (x, t) = — (m 2 — V 2 ,) 4> (x, t) . 



(41) 
(42) 



Formally, the solutions are represented as 

(f)(x,t) = cos yt\J m 2 — (j) (x, 0) + sin \ t\J m 2 — —j= ^ =^ tt (x, 0) 

n(x,t) = 



sin tJm 2 - V 2 J m 2 - V 2 <p (x, 0) + cos [tJm 2 - V 2 ) tt (x, 0) . 



(43) 
(44) 
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From these, it follows by the functional chain rule that 

dx (" <*■ t] sd^j ~ (( m2 " v ') * {x ' ") ^fcy ) (45) 



for any time t. 

Consider the free-field Moyal evolution equation for a generic (not necessarily energy-*-genfunctional) 
WF, corresponding to (|25|) , 



W =-/*(„ W ^ - ♦ (x) ( m > - V«) ^ W. (46) 

The solution is 

W[<j>,v;t] = W[<j>{-t),*{-t);ti\- (47) 

Adapting the method of characteristics for first-order equations to a functional context, one may simply 
check this solution again using the chain rule for functional derivatives, and the field equations evolved 
backwards in time as specified: 

a t W[<t,,ir;t] = d,W[<l>(-t),K{-t);0] 

dx i 9 ' 4, <*■ - 4) wht) + *• {x ' _() ) w l * H) ' ' H) 1 01 

" V & (' (l) 4"('° , " V ')* W *4)) ,i,[ *' l;f| ' (48) 

The quantum Wigner Functional for free fields time-evolves along classical field configurations. In complete 
analogy to the interaction representation for single particle quantum mechanics, (^), the perturbative 
series in the interaction Hamiltonian (written as a ^-function of fields) is then defined in terms of convective 
(time-evolved free field) variables <J?,II: 

ih^W[^,U;t] =H I [$,iq*W[$,Il;t] - W[*,II; *]*#/[*, II| . (49) 

In ref [||, a transformation function T was introduced to accommodate arbitrary canonical transfor- 
mations induced by a generating function F in quantum mechanics, following Dirac. The WF in terms 
of the canonically transformed variables is obtained by convolving with this transformation function. In 
complete analogy, in scalar field theory, cf [J1J], given a canonical transformation from field variables (ft, it 
to variables if, w effected by a generating functional F\<j), cp], one may deduce that the WF in terms of the 
canonically transformed field variables is 



W\6M 



dipdzu 



2tt 



T[4>,ir,<p,w] W[ip,w], (50) 
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where 

T[<t>,ir;ip,m] 



drjdp 



2tt 



expi (f[</> + -ry, cp + -p] - iF*[<f> - -r),ip - -p] + J dx(m (x) p (x) - tt (x) rj (x)) 

(51) 



For example, the generating functional for free field duality between a two-dimensional spacetime scalar 
(p and a pseudoscalar <f> is 

F[cf>,ip} = [dx 4>d x v, (52) 



so it yields the classical canonical transformations 

5 



w = -^F = d x 4>. 

dip 



(53) 



7T = —F = d x cp, 

After some computation, it follows that 

T[<f>, tt; <p, w] = [2tt] 5 [d x ip - tt] 5[w - d x cf>] . (54) 
The ensuing relation between the respective dual Wigner functionals is then quite simple: 



W[<f>,ir] = W 



J 7T, 8 X 



(55) 



to its dual counterpart, [14J, generated by 



A less exceptional example is the canonical transformation from the chiral a-model in two dimensions 

(56) 



F= dx <p l J l (cp) , 



where 



T[p] = Jl-^ <),. • r : ^p J ih/ 



(57) 



is the spatial component of the right (V+A) current. The resulting transformation functional is more 
involved than above, 



[dp] 5 



^-\rU + \p)-\r(p-\p 



X exp ( / / d.v<.>' ( ./' ( p + -p \ - J 1 ftp - -p ) j | <>xp ( / / d.r—p ) . (58) 
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